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A method is derived for calculating the pairing kernel in exchange mediated su- 
perconductors including matrix element effects. Various models for the interaction 
vertex are considered, including spin exchange, orbital exchange, and quadrupo- 
lar exchange. As an example, this formalism is applied to UPt% using relativistic 
wavefunctions from a local density band calculation. 
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Much work has been done over the past decade in trying to understand the microscopic 
basis for superconductivity in heavy fermion metals, and more recently in the copper ox- 
ides. In the heavy fermion case particularly, a large body of experimental data points to the 
presence of a non-trivial order parameter reminiscent of 3 i/e. Because of this, many theo- 
retists working on this subject, including the author, have borrowed techniques successful 
for ^He and applied them to the heavy fermion problem. In particular, a favorite mecha- 
nism being explored is pairing due to spin fluctuations. In the calculations performed, one 
uses some model for the dynamic susceptibility and solves a gap equation. The solution is 
determined by the k dependence of the susceptibility and the Fermi surface 1 ' 2 ' 3 ' 4 . Similar 
calculations have been recently performed on copper oxide materials 5 . 

The deficiency of this approach is obvious. It ignores matrix element effects, that is, the 
k dependence enters only through the energy dispersion and not from the wavefunctions. 
In the analogous magnetic problem, the magnetism would be determined solely by nesting 
effects. Given the commensurate nature of the magnetic order in heavy fermions, nesting 
is an unlikely determinator for magnetism, and thus, one would suspect, for supercon- 
ductivity also. This was realized from the beginning of theoretical work on this subject 6 . 
Despite this, little work has been done along these lines, a noticable exception being a 
paper by Appel and Hertel 7 , where a real space approach to pairing in heavy fermions is 
advocated. 

In this paper, following the spirit of Refs. 6 and 7, the author develops a formalism for 
calculating the pairing kernel in exchange mediated superconductors including relativistic 
matrix elements determined from band structure calculations. Various models for the 
pairing interaction are considered, including spin, orbital, and quadrupolar exchange. As 
an example, this method is applied to the case of UPt^ : the best studied of the heavy 
fermion superconductors. 

The pairing kernel to be evalulated is of the form < k, — k\V\k', — k' > where V is the 
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pairing interaction (note that the band index is implicitly included in the defintion of k). 
This form reduces to those considered in Refs. 1-5 if \k > is taken to be a free electron 
plane wave and V a function of r-r'. In this paper, a different approach is taken by using the 
results of local density band calculations for \k >. Moreover, since it is assumed that the 
physics is primarily determined by the f electrons, an approximation is made of assuming 
that V acts only on the J=5/2 f states on the f atom site. V is further approximated 
by asuming it has just three values corresponding to whether the two electrons are on the 
same site (U), on near neighbor sites (J), and on next near neighbor sites (JN). Finally, the 
dependence on orbital quantum number, fx (where \i ranges from -5/2 to 5/2), is taken into 
account by the exchange process, which is denoted by X, with V taking the form V^ R ,X^Xj 
where i,j represent Cartesian indices (x,y,z), R,R' denote the positions of f atom sites, and 
X will be some appropriate analogue of the Pauli spin matrix, a. Given this, the pairing 
kernel is now of the form 

E V^KklX^k'yuK-klX^-k'y^ (1) 

R,R',i,j 

where < \ >r denotes an integral over a Wigner-Seitz sphere centered at R. To evalulate 
this, note that the band structure wavefunction is of the form 

\k>=^2a^ ktR \/j,> R (2) 

where |/z > is the appropriate combination of spherical harmonics and spinors times the 
f electron radial function. Since X has no dependence on radial coordinate, the radial 
integrals factor out of the problem, except that since an LMTO band method is used, the 
wavefunction has a radial term and its energy derivative as expansion functions. For muta- 
tional convenience, this complication is ignored, but is taken into account when evaluating 
the matrix elements. In a band calculation, \k > is just generated in one irreducible wedge 
of the zone and only on the sites inside the primitive cell. The rest of the information for 
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\k > can be obtained by use of Bloch's theorem plus by application of the group operations 
on \k >. This is somewhat complicated for UPt% since the space group is non-symmorphic 
with two U atoms in the primitive cell. This means that some group operations must be 
followed by non-primitive translations and others lead to the interchange of the two U 
sites. The effect of group operations on \ji > is well known, for UPt% they transform as 
i~7, .Tg, and i~g. This has been tabulated in Ref. 7, for instance, but the author rederived 
all relations used in the programs as a safety check. Finally, in the approximation used in 
this paper, Bloch's theorem will appear by factors of the form in the matrix elements 
where R is the site position of a U ion. These cancel one another for the on-site case, but 
must be carefully kept track of for near neighbor and next near neighbor terms. 

The nature of the operator X will now be discussed. In the ^He problem, it is taken 
to be the Pauli spin matrix. The simplest generalization for the current case is to replace 
a by the pseudo-spin matrix r, where r has the same effect on k, PTk (or Pk, Tk) that a 
has on up and down spinors. Note that P is the parity operator, and T the time reversal 
operator. Further generalizations can be made by replacing r by J, the total angular 
momentum operator. Further, if one is interested in quadrupolar interactions, thought by 
some to be of fundamental significance in the heavy fermion problem^, then one can in 
turn replace J by 0\ = [3jf - J(J + 1)]/V5, 0\ = [J Z (J+ + J-) + (J+ + J-)J z ]/2, or 
0\ = \J\ + J 2 ]/2. 

To determine the appropriate symmetry of the gap function, one must project the pair- 
ing kernel onto each group representation. For UPt% these are i~i through Ig, equivalent to 
the more common "chemical" notation of A\, A% B\, E>2, E2, E\. Moreover, a projection 
must also be done on odd and even parity. For even parity, the pair function notation 
\k,—k > is replaced by \k,Tk > —\PTk,Pk >. The odd parity case involves a three 
component vector. An assumption will be made in this case. Experimentally, it appears 
from neutron scattering data 9 that the moments in UPt^ are confined to the basal plane. 
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If a similar effect occurs for the Cooper pair moments, then one expects only the n d z n 
component of the odd parity order parameter vector to survive. Such an order parameter 
is consistent with a recent analysis of H c 2 anisotropy in UPt^ 10 . This projection automat- 
ically occurs in the r model if one assumes that V xx = V yy with the rest of the V lJ equal 
to zero. For the other models, mixing will occur. It is assumed, though, that there is some 
term in the pairing Hamiltonian which also acts in these cases to project onto d z . For the 
relativistic case, d z is \k,Tk > + \PTk,Pk >. 

To summarize, the following pair intractions will be explored: T x T x + TyTy, JxJx + JyJy, 
0§0§, 0\0\, 0\0\. \k,-k> will be taken to be \k,Tk > T\PTk, Pk > where the upper 
sign is for even parity and the lower one for d z odd parity. The weak-coupling gap equation 
to be solved is 

4fe = A -1 ZXfc'Afe' (3) 
k' 

where V^ y is the pairing kernel discussed above and A -1 , the inverse coupling constant, is 
ln(1.13u> c /T c ) with T c the transition temperature and u c some cut-off energy. The sum is 
carried out over all k vectors on the Fermi surface. For UPt^ : a 137 k vector mesh is used 
on the Fermi surface in the irreducible wedge, noting that there are 24 group operations in 
this case. The k vectors are weighted by the appropriate denstiy of states factors derived 
from a tetrahedron decomposition of the Brillouin zone. 

In Table 1, the results of the above formalism for UPt% are summarized. In column 2, 
an on-site interaction of unit value is looked at. The sign in each case is chosen opposite 
to that which would yield a nodeless A\ g (s-wave) solution. This sign is positive for the r, 
J, and 0\ cases, but is negative for the 0\ and 0\ cases. Despite the "repulsive" value 
of the sign, solutions are indeed found, due to the momentum dependence of the matrix 
elements. The solution with the largest coupling constant either has A\ g symmetry or is 
odd parity (for the r case, the interaction is repulsive for all even parity states). Addition 
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of a repulsive momentum independent constant to the gap equation did not suppress the 
A\g solution since this solution already has a complicated nodal structure. In fact, the 
stability of the A\ g solution over all other even parity solutions is most likely due to this 
complicated nodal structure, since the nodal structure of the other representations is fixed 
by symmetry and therefore not free to adjust itself to minimize the free energy. 

In the final four columns, results for ferromagnetic (positive) and antiferromagnetic 
(negative) unit values of J (near neighbor) and JN (next near neighbor) interactions are 
tabulated. Again, the solution with the largest coupling constant either has A\ g symmetry 
or is odd parity. An interesting point to remark on is that for the r and J cases, one finds 
that the ferromagnetic sign supports odd parity solutions whereas the antiferromagnetic 
sign supports even parity (for the r case, the ferro case is repulsive for all even parity states 
and the antiferro case is repulsive for all odd parity states). This is in accord with earlier 
wisdom on this subject 1 but did not occur for realistic calculations on U Pt$ which ignored 
matrix element effects 2 . This illustrates the crucial importance of matrix elements when 
evalulating the pairing kernel. Another interesting point is that quadrupolar interactions 
prefer A\ g , regardless of the sign, in ten of twelve cases. 

The most favored explanation for the large body of experimental data on UPt% is 
that the order parameter is from a two-dimensional group representation with line nodes 
orientated perpendicular to the c axis 11 . In the current context, this could be realized 
by E\ g or E^ui the latter being consistent with the H c 2 anisotropy analysis 10 . Neither 
state is favored by any of the pair interactions analyzed in the current paper (the same 
was found for earlier work which ignored matrix element effects 2 ). The most favored two 
dimensional rep in the current work is E\ u . This solution has point nodes. On the other 
hand, arguments have been put forth that the observed line node structure in U Pt% may 
not be due to the symmetry of the order parameter, but rather to gaplessness due to the 
normal state self-energy 12 . More work will certainly be necessary before the actual form 
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of the order parameter is unambiguously determined. 

An alternate explanation of the data would be a near degeneracy of two different 
representations 13 . This picture is certainly supported by the current work, which indicates 
several such cases which involve one single dimensional rep and one two dimensional rep. 
The picture advocated by Machida and co-workers 14 which assumes a near degeneracy of 
the three components of an odd parity solution is not consistent with the treatment of this 
paper where strong spin-orbit effects prohibit such a degeneracy. 

The author concludes with some remarks on strong coupling effects. The author has 
performed calculations on UPt% including the full momentum and frequency dependence 
of the pairing interaction and treating all k vectors and f bands (instead of restricting to 
the Fermi surface) as recently advocated by Monthoux and Pines 5 , but ignoring matrix 
element effects. Although this treatment leads to a renormalization of the value of T c , no 
change in the ordering of the solutions with respect to coupling constant was seen. For 
that reason, one would suspect that strong coupling effects would not alter the ordering of 
solutions seen here, although it is conceivable that variation of the matrix elements when 
going off the Fermi surface could alter the results. The author plans to test such effects 
in the future if the effort seems warranted. An alternate strong coupling effect to consider 
is feedback to the pair potential caused by changes induced when the metal enters the 
superconducting phase. These effects are crucial in the case of 3 ife, but for the heavy 
fermion case, there has been little evidence for a change in the dynamic susceptibility 
below T c . Until such changes are identified, it would be premature to invoke such effects 
in the current formalism. 

As a final thought, the author would like to emphasize the generality of the above 
method. Currently, this has been applied to a phenomenological and very simple interac- 
tion vertex for UPt% as an example. Once a more microscopic theory is available, then the 
resulting calculation for properly testing such a theory must follow a similar formalism as 
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expounded in this paper. This statement will most likely apply to other cases which might 
involve exchange mediated interactions, such as the copper oxide materials. 

In summary, a new method is derived for evaluating the pairing kernel in exchange me- 
diated superconductors utilizing matrix elements determined from relativistic band struc- 
ture calculations. As an example, several models for the pairing interaction for UPt% 
were considered, yielding a rich variety of solutions. Further work, both experimental and 
theoretical, will be necessary before it can be determined what relevance the determined 
solutions have to the experimental situation. 

This work was supported by U.S. Department of Energy, Office of Basic Energy Sci- 
ences, under Contract No. W-31-109-ENG-38. 
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Table 1. Pairing coupling constants for various group representations for UPt%. U repre- 
sents an on-site interaction of unit value, with the sign chosen opposite that which would 
give a nodeless s-wave solution. N and NN represent near neighbor and next near neigh- 
bor interactions of unit value with either a positive sign (F) or negative sign (AF). The 
interaction vertex is denoted by r for t x t x + T y T y , J for J x Jx + JyJy, 2 for 2 2 , 0\ for 
0\0\, and 0\ for 0\0\. Listed are the group representations with the largest coupling 
constants (with the coupling constant in parenthesis) . 



Case 


u 


F-N 


AF-N 


F-NN 


AF-NN 


r 


Sl„(0.015) 


£l u (0.040) 


A lg (0.064) 


S lu (0.047) 


A lg (o.m) 




A 2u (0.009) 


Si u (0.036) 


E 2g (0M7) 


^1 U (0.046) 


E 2g (0.070) 


J 


A lff (0.058) 


£?l„(0.188) 


A lg (0.209) 


S lu (0.280) 


A lg (0.728) 




Sl u (0.047) 


Sl«(0.169) 


£ 2g (0.172) 


E iu (0.236) 
A lff (0.207) 


E 2g (0.397) 




A 2m (0.247) 


Ai g (0.838) 


A lg (1.33 ) 


A lg (2.32 ) 


A 2?i (0.484) 




£?l„(0.171) 


A 2u (0.699) 


E 2g (0.600) 


£ lu (0.635) 


S lu (0.334) 




£ 2ff (0.161) 


£ lu (0.392) 


E lu (0A13) 


S lu (0.629) 


A l5 (0.328) 
^ln(0-323) 


o\ 


A lff (0.157) 


A lg (0AS8) 


A lg (0.361) 


A lg (1.60 ) 


A lg (0.547) 




Bi u (0.121) 


A 2u (0.216) 


E 2g (0.227) 


£ 2ff (0.486) 


E lu (0.323) 




£ 2ff (0.102) 


Sl„(0.215) 
£?l„(0.194) 


A 2u (0.186) 




S lu (0.272) 


ol 


A lg (0.U5) 


A lg (0.374) 


Ai ff (0.236) 


A lff (0.381) 


A lg (1.05 ) 




Bi„(0.141) 


El 9 (0.100) 


El u (0.235) 
E 2g (0A66) 


E 25 (0.174) 
Sl„(0.163) 


£ 2g (0.571) 
Sl u (0.392) 



